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ABSTRACT 

Observations of outflows associated with pre- main- sequence stars reveal details about morphol- 
ogy, binarity and evolutionary states of young stellar objects. We present molecular line data from 
the Berkeley-Illinois-Maryland Association array and Five Colleges Radio Astronomical Observatory 
toward the regions containing the Herbig Ae/Be stars LkHa 198 and LkHa 225S. Single dish obser- 
vations of 12 CO J = 1 - 0, 13 CO J = 1 - 0, N 2 H+ J = 1 - and CS J = 2 - 1 were made over 
a field of 4.3' x4.3' for each species. 12 CO J = 1 — data from FCRAO were combined with high 
resolution BIMA array data to achieve a naturally- weighted synthesized beam of 6.75" x 5.5" toward 
LkHa 198 and 5.7" x 3.95" toward LkHa 225S, representing resolution improvements of factors of 
approximately 10 and 5 over existing data. By using uniform weighting, we achieved another factor of 
two improvement. The outflow around LkHa 198 resolves into at least four outflows, none of which are 
centered on LkHa 198-IR, but even at our resolution, we cannot exclude the possibility of an outflow 
associated with this source. In the LkHa 225S region, we find evidence for two outflows associated 
with LkHa 225S itself and a third outflow is likely driven by this source. Identification of the driving 
sources is still resolution-limited and is also complicated by the presence of three clouds along the line 
of sight toward the Cygnus molecular cloud. 13 CO J — 1 — is present in the environments of both 
stars along with cold, dense gas as traced by CS J = 2 — 1 and (in LkHa 225S) N 2 H + J = 1 — 0. 
No 2.6 mm continuum is detected in either region in relatively shallow maps compared to existing 
continuum observations. 

Subject headings: ISM: clouds, molecules — ISM: individual (LkHa 198) — ISM: individual (LkHa 
225S) — stars: formation — millimeter 



1. INTRODUCTION 

In addition to evidence for mass infall (accretion), 
many pre-main-sequence stars also possess prominent 
outflows and jets. These jets are believed to play a cru- 
cial role in removing angular momentum from circumstel- 
lar d isks, allowing accretion to continue ([Stahler fc Pallal 
12004( 1. Outflows can extend many parsecs in length, 
and can serve as "fossil records" of the accretion pro- 
cess, providing indications of temporal variability of ac- 
cretion or_^rece^sion_oXoutflow axes over thousands of 
years ([McGroartv fc Ravll2004l ). Furthermore, outflows 
inject turbulence into molecular clouds and may regulate 
overall star formation effic iency (|Matzer fc McKeell2000l : 
IChrvsostomou et al.ll2005[ ). 

Outflows from young stars may be observed 
both through shock-excited e mission lines such 
as [SHI 6716 6731 A (e.g., iFinkenzellerl 119851 : 
ICorcoran fc Ravi Il998h and [Felll 1.257, 1.644 urn 
(jGraham. Wright fc Longmord 119871 : iPvo et all I20Q3I ). 
and also through millimeter rotational lines of CO 
and other molecules (e.g., ISchmid-Burgk et al.l 1990; 
lArce fc Sargentl 120061 ). Because star formation is a 
complex and turbulent process, high angular resolution 
observations (i.e., submillimeter or millimeter interfer- 
ometry) are usually required in order to disentangle 
the various motions of molecular gas around young 
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stars (infall, outflow, rotation, or turbulence). Such 
observations have been used to probe the kinematics of 
outflows from a wide range of YSOs, i ncluding T Tauri 
stars such as HH 30 ( Pet v et all T2006) and also deeply 
embedded C l ass I protostars such as in AFGL 5142 
(|Zhang et al.l l2007t ), which is a protocluster displaying 
at least three molecular outflows originating from at 
least half a dozen embedded massive protostars. 

In this paper we consider molecular outflows in an in- 
termediate case between these two extremes, by investi- 
gating outflows originating from groups of Herbig Ae/Be 
(HAeBe) stars. The HAeBe stars are pre-main-sequence 
stars of intermediate-mass, i.e., 1.5 < M/Mq < 10 
(jHerbid [l960t I Waters fc Waelkinsl [l998h . which exhibit 
infrared and millimeter excesses indicative of circum- 
stellar dust. While the majority of outflows studied to 
date originate from T Tauri stars, many Herbig Ae/Be 
stars are known to launch outflows as well (|Mundt fc Ravi 
fl99l. in some c ases extending several parsecs in length 
(|McGroartv et all I2004T) . HAeBe stars are often found 
m groups (jTesti et all I1998L 11997( 1. making high angu- 
lar resolution observations critical for resolving the out- 
flows associated with individual stars. These intermedi- 
ate mass stars are also rarer than later type T Tauri stars 
due to the steepness of the Initial Mass Function, which 
means they are, on average, found at greater distances 
than their lower mass counterparts. Therefore, higher 
angular resolutions are required to resolve their outflows 
compared to closer, lower mass T Tauris. 

We have searched for wide field and high resolution 
emission from 12 CO J = 1—0 in the environments around 
the HAeBe stars LkHa 198-IR and LkRa 225S based on 
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TABLE 1 
FCRAO Observational Summary 



Source 


Species/Transition 


Frequency 


Antenna 


Resolution 


Sensitivity per Channel 






[GHz] 


Efficiency 


[km s- 1 ] 


[K] 


LkHa 198 


12 CO J = 1-0 


115.21720 


0.45 


0.06349 


0.35 




13 CO J = 1-0 


110.20135 


0.49 


0.06642 


0.15 




N 2 H+ J = 1-0 


93.17378 


0.52 


0.07855 


0.09 




CS J = 2 - 1 


97.98095 


0.52 


0.07470 


0.09 


LkHa 225S 


^CO J = 1-0 


115.21720 


0.45 


0.06350 


0.40 




13 CO J = 1-0 


110.20135 


0.49 


0.06642 


0.16 




N 2 H+ J = 1 — 


93.17378 


0.52 


0.07470 


0.09 




CS J = 2 - 1 


97.98095 


0.52 


0.07856 


0.09 



intri guing; results from near- IR AO imaging polarimetry 
by iPerrin et al.1 (|2004af ). In addition, we have obtained 



single dish spectra of 13 CO J = 1 — 0, CS J = 2 — 1 and 
N2H + J = 1 — toward both sources. We present here 
millimeter observations of the LkHa 198 and LkHa 225S 
stars and their nearby group members. Our high-angular 
resolution observations resolve multiple outflows in both 
of these regions. The observations and data reduction 
techniques are described in §[2j Results for LkHa 198 and 
LkHa 225S and their surrounding regions are presented 
in §[3] and §21 respectively. Our findings are summarized 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. FCRAO Data 

We obtained 12 CO J = 1 - 0, 13 CO J = 1 - 0, 
CS J = 2 - 1 and N 2 H+ J = 1 - data from the 
14 m Five Colleges Radio Astronomical Observatory 
(FCRAO) using the array SEQUOIA on 2004 Novem- 
ber 10. The telescope was pointed toward LkHa 198 at 
«J2000 = 00 h ll m 25*97, (5j2ooo = +59°49 / 29 / /l and LkHa 
225S at a J2 ooo = 20 h 20 m 30*65, £ J20 oo = +41°21 / 25 / /50. 
A fully-sampled image covering a field of view 4.3' x 4.3' 
was made of four spectral line transitions. The tran- 
sitions 12 CO J = 1 - and 13 CO J = 1 - were ob- 
served simultaneously; similarly, CS J = 2 — 1 and N 2 H + 
J = 1 — were obtained in a single observation. A sum- 
mary of the FCRAO observations is shown in Table [TJ 
Each observation was 24 minutes in duration. Calibrated 
T^ data are delivered from the telescope. We converted 
these data to main beam temperatures by division of the 
appropriate antenna efficiency (see Table [1]). The conver- 
sion to Jy beam -1 was done by deriving the conversion 
factor based on the observing wavelength and telescope 
beam: 



Jy 

K 



13.6 



A [mm]' 



mm max 



(i) 



where 6^ in and Q max are the minor and major axes of 
the telescope beam. For single dish millimeter telescopes, 
the axes are equal and equivalent to the diffraction limit 
of the antenna. For FCRAO, the conversion factor is 
23.6 Jy/K. The beamsize (resolution) is 46.7" for 12 CO 
J = 1-0, 48.9" for 13 CO J = 1-0, 55.0" for CS 
J = 2 - 1, and 57.5" for N 2 H+ J = 1 - 0. The data 
were taken in the 25 MHz mode. The rms sensitivity per 
channel is shown in Table [1] for each transition. 

2.2. BIMA Interferometric Data 



We observed the regions surrounding the two Her- 
big Ae/Be multiple systems in single pointings with the 
Berkeley -Illinois-Maryland A ssociation (BIMA) interfer- 
ometer 3 ([Welch et al. Ill 9961 ) in Hat Creek, CA. The data 
were taken in 2003 October (1 track each) followed by 
multiple observations over the period of 2004 March to 
2004 June (see Table |2j). For both regions, we used the 
same target coordinates as for the FCRAO observations. 

Three configurations of the ten 6.1 m antennas were 
used to observe the 12 CO J = 1 — line and continuum 
at 2.7 mm. An A- array track and a D-array track toward 
LkHa 225 were not used due to poor phase correlations. 
Two B- array tracks toward LkHa 198 were rejected for 
the same reason. The D-array, C-array and B-array data 
have projected baselines between 2-13 kA, 2-33 kA and 
3-74 kA, respectively. Table [2] contains the sensitivities 
achieved per useable track. 

The 12 CO J = 1 — line was observed utilizing the 
digital correlator to record the line in four bands, config- 
ured in redundant bandwidths of 12.5 MHz (256 chan- 
nels) and 50 MHz (128 channels), giving spectral res- 
olutions of 0.127 km s _1 and 1.016 km s _1 for LkHa 
198. The same configuration was used for the C-array 
track of LkHa 225S. For the 2004 tracks toward LkHa 
225S, we allocated 256 channels each across bandwidths 
of 12.5 MHz and 25 MHz, resulting in velocity resolu- 
tions of 0.127 km s _1 and 0.254 km s _1 , respectively. 
LkHa 198 was also observed in some LkHa 225S tracks, 
when the latter source was too high in elevation to be 
observed. The lower sideband was used to measure the 
continuum emission at 2.6 mm in bandwidths (125 and 
75 MHz for the respective configurations) limited by the 
spectral line correlator configuration. The typical BIMA 
array continuum bandwidth was 800 MHz. 

Phase and amplitude variations were calibrated by ob- 
serving quasars at small angular separations from the 
source fields: 2015+372 and 0102+584 (for LkHa 225S 
and LkHa 198 respectively) approximately every 30 min- 
utes. The adopted fluxes of these quasars were epoch de- 
pendent and measured against observations of the planet 
Uranus or the quasar 3C454.3 when possible. The cal- 
ibration was performed using the MIRIAD (Multichan- 
nel Image Reconstruction, Ima ge Analysis and Display; 
LSault. Teuben fc Wright] (jl995h ) task MSELFCAL. Ab- 
solute flux calibration was done using Uranus when ob- 
served or by using the derived fluxes of gain calibrators 

3 The BIMA Array was operated by the Berkeley-Illinois- 
Maryland Association under funding from the National Science 
Foundation. 
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TABLE 2 

Summary of 12 CO J =1-0 Observations at 115.21720 GHz 



Source 


Date 


Array 


Resolution 


Sensitivity per Channel 






Configuration 


[km s J 


[Jy beam -1 ] 


LkHa 198 


1 March 2004 


B 


0.127 


4.6 




4 March 2004 


B 


0.127 


5.7 




18 March 2004 


B 


0.127 


6.6 




29 March 2004 


B 


0.127 


3.4 




26 April 2004 


C 


0.127 


1.9 




1 June 2004 


D 


0.127 


1.8 


LkHa 225S 


28 October 2003 


C 


0.127 


2.7 




4 March 2004 


B 


0.254 


0.85 




21 March 2004 


B 


0.254 


1.6 




31 May 2004 


D 


0.127 


5.9 



TABLE 3 

Final Parameters of Combined FCRAO + BIMA 12 CO J = 1 - Maps 



Source 


Channel Width 
[km s- 1 ] 


Final Beam 




Weighting 


rms per channel 
[Jy beam -1 ] 


LkHa 198 


0.127 


6.75" x 5.47" at 


-87° 


natural 


0.66 




0.127 


3.39" x 2.37" at 


-82.5° 


uniform 


1.4 


LkHa 225S 


0.254 


5.67" x 3.95" at 


-79° 


natural 


0.45 




0.254 


2.87" x 2.29" at 


-89° 


uniform 


0.70 



during the same epoch as our observations (utilizing the 
catalogue of fluxes at "plot_swflux" on the BIMA web- 
site 4 ). 

Subsequent processing of the data, including the com- 
bination of data from different configurations, was done 
with MIRIAD. Images were produced using MIRIAD's 
INVERT algorithm with natural weighting and then 
by variation of the "robust" parameter to produce a 
uniform-weighted images optimized for both signal-to- 
noise and spatial resolution. The images were then de- 
convolved using CLEAN and then convolved with the 
clean beam using RESTOR. The resulting resolutions 
and sensitivities of the BIMA data are shown in Table 
[2j The resulting noise levels in naturally-weighted im- 
ages are 0.9 Jy beam -1 in 0.127 km s -1 channels for the 
12 CO J = 1 — line emission from LkHa 198 and 0.75 Jy 
beam -1 in 0.254 km s _1 channels for LkHa 225S. The 
rms levels in the continuum maps are 18.5 and 17.0 mJy 
beam -1 toward LkHa 198 and LkHa 225S respectively. 
We achieved naturally-weighted FWHM spatial resolu- 
tions of 7.4" x 6.0" and 5.5" x 4.1" for the continuum 
maps of the two fields. 

2.3. Combination of FCRAO and BIMA data 

The 12 CO J = 1 - data from FCRAO and the BIMA 
array were combined to produce total power data cubes. 
The 12 CO J = 1 — emission is strongly detected toward 
both sources, so we have used the highest possible spec- 
tral resolution data for combination with the FCRAO 
data. In the observation of LkHa 225S, one spectral win- 
dow exhibited evidence of corrupt spectra; therefore, we 
utilized a redundant observation in a second spectral win- 
dow for that source. The FCRAO data were reordered 
to match the BIMA axes (x,y,v) and rescaled to Janskys. 
The MIRIAD task REGRID was then used to interpolate 

4 http://bima.astro.umd.edu 



the FCRAO data onto the BIMA d ata cube. We used 
the li near combination algorithm of IStanimirovic et al.l 
(1999) to combine the naturally- weighted BIMA array 
dirty map with the single dish map, using the ratio of 
the beams as the weighting factor. The combined image 
was then deconvolved and convolved with the clean beam 
to reconstruct the total 12 CO J = 1 — emission. Maps 
were created with natural and uniform weighting. The 
resultant resolutions and sensitivities are summarized in 
Table 03 The primary beam of the BIMA array anten- 
nas at 115 GHz is 1.8', which is the area of the resultant 
combined data maps. 

3. LkHa 198 AND ENVIRONS 

LkHa 198 (V* V633 Cas, IRAS 00087+5833) is an 
HAeBe star located at ~ 600 pc based on the inter - 
stellar extinction l aw dRacind 1 19681 : IChavarria-K.lll985h . 
iHillenbrand et all (|l992l ) categ orize LkRa 198 as an A5 
pre-main sequence star while iHernandez et al.l (2004) 
find a spectral type of B9 ±2.5. LkHa 198 is classi- 
fied as a Group I I HAeBe star in the nomenclature of 
Hill enbrand et all Jl992), 

meaning: that the SED rises 
with wavelength in the infrared. Group II HAeBe stars 
have significant circumstellar envelopes which reprocess 
stellar radiation at long wavelengths. A strong CO bipo- 
lar outflow is associated with this region, but it is not 
centered on LkHa 1 98, but rather to the northwest 
(jCanto et al.lll984l ). iStrom et all (|l986f ) detected evi- 
dence of a blue-shifted optical flow extending to the east- 
southeast of LkHa 198, which did not coincide with the 
red-shifted CO emission in that direction. An infrared 
companion (LkHa 198-I R) was identified a pproximately 
6" north of LkRa 198 by lLagage etaD (|1993[ ). A submil- 
limeter source (LkHa 198-mm) wa s identified 19" north- 
west o f LkHa 198 at 800 /xm by ISandell fc Weintraubl 
(|l994l ) who also found evidence for CO outflows in the 
three sources LkHa 198, LkHa 198 mm and V376 Cas. 
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Fig. 1. — Left: The greyscale shows the zeroth moment over the entire range of CO emission (-4 to 4 km s -1 ) of the combined FCRAO 
+ BIMA data toward the region surrounding LkHo; 198 (in units of Jy beam -1 km s _1 ). Optically visible HAeBe stars are marked as 
stars, while the embedded millimeter source is indicated by the square. Contours illustrate the emission to the blue (—4 to —1 km s _1 ) and 
to the red (1 to 4 km s _1 ) of the peak. The contour levels range from 5.2 to 26 Jy beam -1 km s _1 in steps of 2.6 Jy beam -1 km s -1 . The 
resolution is 6.75" x 5. 5". At 600 pc, the beam corresponds to a linear scale of 4000 x 3300 AU. The box indicates the field of view of the 
uniform- weighted image at right. Right: A high resolution datacube is produced from uniformly- weighted BIMA array data combined with 
FCRAO data. The greyscale shows zeroth moment over the entire range of CO emission (—4.0 to 4.2 km s -1 ). The blue and red contours 
show the zeroth moments over the ranges from —4.0 to —2.1 km s -1 and 2.3 to 4.2 km s -1 , respectively. The contour levels range from 
1.8 to 4.2 Jy beam -1 km s -1 in steps of 0.6 Jy beam -1 km s _1 . The resolution of 3.4" x 2.4" (2040 x 1440 AU) highlights primarily 
small features in the outflowing gas. This map shows that no significant r ed- or blue-sh i fted em ission is associated with the LkHo; 198-IR 
source. For comparison, the polarized intensity in H-band as measured by Perrin et a l. (2 004a ) is shown in the inset. The LkHo; 198 and 
LkHa 198-IR sources are marked on this image. 



This source is another HAeBe star in the region, located 
about 40" north of LkHa 198. LkHa 198 mm was also 
str ongly detected i n a 1. 3 mm observation of the region 
by iHenning et al.l (|l998l ) . High resolution interferome- 
try toward LkHa 198 and the surrounding region at 2.7 
mm failed to detect Lk Ha 198, but LkHo _ 198-m m was 
detected and resolved (|Di Francesco et aDll997f ). High 
resolution la ser guide sta r adapt ive optics imaging po- 
larimetry by iPerrin et al.l (|2QQ4af ) revealed a centrosym- 
metric polarimetric pattern consistent with a conical cav- 
ity. They also detected a "cometary" structure in JHK 
intensity images. Their data suggest that this structure 
is associated with LkHa 198-IR. 

3.1. Observational Results 

3.1.1. Moment and Channel Maps 

Figure [1] shows the zeroth moment of the combined 
map of 12 CO J = 1 — emission toward the LkHa 198 
region from the combination of channels from —4 to 4 km 
s _1 and the distribution of emission at velocities blue- 
ward and red ward of the systemi c LSR source velocity, 
—0.7 km s _1 ([Canto et al.l 1 1981 ). The moment maps 
were created by using only channels with flux densities 
exceeding 1 Jy beam -1 , approximately 1.5 times the rms 
level per channel. The velocity ranges over which they 
are taken are specified in Figure [TJ The resolution in this 



map is an order of magnitude improved over the 12 CO 
J = 1 — maps of iCanto et all (1984) but covers an 
area limited by the primary beam of the BIMA array 
antennas. Strong 12 CO J = 1 — emission is concen- 
trated around LkHa 198 rather than LkHa 198-IR to 
the north. We have not detected any structures simi- 
lar to the dust extension f rom LkHa 198-IR o bserved in 
near-IR polarimetry data (Per rin et al.ll2QQ4ah . 

The second panel of Figured] shows a higher resolution 
image of the region surrounding LkHa 198 and LkHa 
198-IR produced by uniformly weighting the BIMA ar- 
ray data before combination with the FCRAO map. The 
CO features in this map are concentrated around LkHa 
198, with no significant emission detected around LkHa 
198-IR. Blue- and red-shifted emission peaks are paired 
on either side of LkHa 198. If these each represent indi- 
vidual outflows driven by independent sources, then the 
driving sources are not evident as no sources of strong 
dust emission are located at the central positions be- 
tween the red and blue peaks. Instead, it is more likely 
that the peaks represent parts of an outflow from LkHa 
198 which may be precessing over time. Such precession 
can occur when the outfl ow's driving source is a mem- 
ber of a multiple system. Sm ith et al.l (|2QQ5l ) show that 
LkHa 198 is a spectroscopic binary with a separation of 
60 mas. 
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Fig. 2. — Binned intensity images of 12 CO J = 1 — from combined FCRAO and BIMA array data toward the LkHa 198 region is 
shown in greyscale and contours over 1.016 km s _1 bins. The mean velocity per channel is noted. The greyscale range is common to all 
channels and is shown at right (Jy beam -1 ). The contours are plotted at 20 to 200 a in steps of 20 a where a = 0.64 Jy beam -1 is the 
rms level in a single channel before binning. Sources are marked by asterisks and the square as in Figure ^ 



Figure [2] shows binned channel maps with 1 km s _1 
steps around the V(LSR) velocity of —0.7 km s _1 . In the 
blue-shifted channels, the dominant feature is the arc of 
emission that traces out the edges of the s outheast ellip- 
tical c avity detected in the I-band image oflLagage et all 
()l993h and the J-band image of iKoresko et all (|1997). 
This feature appears more strongly focussed on LkHa 
198 than the adjacent LkHa 198-IR. A secondary blue- 
shifted feature to the north (extending off the map area) 
may be a component of the outflow originating at LkHa 
198-mm. The channels around the cloud's velocity are 
strongly dominated by cloud emission, but the peak as- 
sociated with LkHa 198 persists over most channels. By 
the 1.2 km s _1 channel, a red-shifted outflow component 
can be associated with LkHa 198-mm. The red-shifted 
counterpart to the LkHa 198 outflow appears along the 
same axis. At the highest velocity shown, 3.3 km s _1 , 
components of the two outflows can be distinguished. In- 
terestingly, red-shifted emission to the south and south- 
ea st also seems t o trace out the elliptical cavity detected 
bv lLagage et al.l (|l993l ) in the optical. Interestingly, we 
find no strong 12 CO J = 1 — features associated with 
LkHa 198-IR. 

We note that Figure [T] may seem ini t ially at odds with 
the results of Figure 3 of I Canto et al.l (|l984l ) who report 
single dish detection of 12 CO J = 1 — in a 66" beam. 
The red lobe peak was found to extend to the southeast 
and to be centered on LkH a 198 with a 1.5' off set to the 
blue peak to the northwest. I Canto et al.l ([1 984) did note 
that their observations did not depict lobes symmetric 
about the star, as is more typical for a bipolar outflow, 
but they did interpret the blue and red shifted emission 
as an outflow driven by LkHa 198. It is noteworthy as 



well that blue-shifted optical emission was detec ted to 
the southeast of LkHa 198 bv IStrom etafl (|1986f ). Our 
single-dish FCRA O data have high er angular resolution 
(46.7") than the lCanto etafl (l984) map. We have there- 
fore convolved our FCRAO data to 66" resolution (see 
Figure [3|). When compared to the existing dataset, we 
find a similar morphology in the zeroth moment map, 
with a single prominent red lobe to the southeast and a 
blue lobe to the northwest. The FCRAO plot in Figure 
[3] has been generated using the same range in velocity 
as Figure 3 of Canto et al. and plotted with comparable 
flux density levels. At a resolution of 66", the emission in 
the region is strongly dominated by red-shifted emission 
which is centered on the LkHa 198/198-IR sources. The 
blue-shifted emission i n the map lies appr oximately V to 
the north (as found by ICanto et al.ll 19841 ). The fact that 
the outflow has been acknowledged to be very asymmet- 
ric about the att ributed driving source, either LkHa 19 8 
or LkHa 198-IR ([Canto et al.lll984l : iNakano et al.lll99QT ). 
could point to multiple driving sources, but could also 
be explained if the outflow is not much inclined from the 
line of sight. 

The ambiguity regarding the asymmetric nature of the 
outflow is resolved by the zeroth moments of the BIMA 
array data alone, which are also shown in Figure 03 This 
map illustrates sharply the effect of a ten-fold improve- 
ment in resolution on the interpretation of the number 
and orientation of outflows in a region containing multi- 
ple sources. Even though we are only sensitive to struc- 
tures on limited spatial scales with these data, the fea- 
tures in red and blue emission can be more clearly as- 
sociated with individual objects than is possible based 
on single dish observations. Figure [3] shows that the red 
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Fig. 3. Zeroth moment maps of the LkRa 198 re gion in 12 CO J = 1 - from FCRAO (top left), the BIMA array (bottom left) 
and the combined dataset (bottom right). Figure 3 of I Canto et al.l ([19841 ") is reproduced at the top right. The grey box shows the area 
covered by our maps. The contours of the FCRAO image are produced by integration over the same velocity channels as Canto et al. (i.e., 
—4.7 to —1.5 km s _1 for the blueshifted emission and 0.7 to 3.9 km s _1 for the redshifted emission) and plotted at comparable levels (312 
Jy beam -1 to 468 Jy beam -1 in steps of 52 Jy beam -1 ). Channels with flux densities less than 35 Jy beam -1 were not included. The 
elongated red lobe dominates the emission. The blue lobe of I Canto et al.l (1 19841 ) is just detected at the NW edge of our map, which is less 
extensive in that dir ection. The r e soluti on of the FCRAO data is 46.7" (grey circle), which we have convolved to 66 /x (white circle) to equal 
the resolution of the Canto et al. (1984) map. The BIMA array data are integrated over the same velocity ranges, with only channels with 
flux densities exceeding 1 Jy beam -1 included. The circle denotes the primary beam of BIMA at 115 GHz. This is effectively the edge of 
sensitivity in the BIMA and combined maps. The combined FCRAO+BIMA map is made with the same contour levels as shown in Figure 
^and integrated over the same velocities. The BIMA array and combined maps have a synthesized beam of 6.75 /x x 5.5 /x as indicated at 
the bottom right of those maps. The sources are indicated by asterisks and the square as for Figure [T] 



lobe resolves into distinct features, likely associated with 
LkHa 198 and/or LkHa 198-IR, a distinct red peak as- 
sociated with LkHa 198-mm, and larger scale blue- and 
red-shifted emission in the vicinity of V376 Cas. Whether 
the blue-shifted emission to the southeast is driven by 
LkHa 198 or LkHa 198-IR cannot be easily distinguished 
in this map, since their small separation of 5" is compa- 
rable to our resolution. However, the channel maps of 
Figure [2] do favour LkHa 198 as the driving source for 
both outflows originating in that small area. The im- 
proved resolution of our data explains the coincidence of 
th e blue-shifted optic al flow defined by the [SII] emission 
of IStrom et"afl (I1986D within t he red -shifted CO outflow 
lobe detected bv ICanto et all (|l984t ). 

With the higher resolution data in hand, we can in- 
terpret the single dish data as follows. The red-shifted 
emission originating from the LkHa 198/ 198-IR objects 
completely dominates the emission near the four sources, 
with the weaker blue lobe shifted north because it orig- 
inates from one or both of V376 Cas or LkHa 198- 
mm. This general 12 CO J = 1 — morphology in Fig- 
ur e [T] is discernable from the CO J — 2 — 1 emission 
of ISandell fc Weintraubl (|l994f ) observed with 20" res- 
olution with the JCMT, which shows that blue-shifted 
and red-shifted emission overlap southeast of the LkHa 
198/198-IR sources and that secondary blue and red- 



shifted peaks are seen > V north of these objects. How- 
ever, with the superior resolution of the BIMA array (and 
the removal of large scale emission from the maps), the 
outflow structure can be resolved. 

3.1.2. Spectra 

In addition to the images presented above, we can con- 
struct spectra toward the known point sources in the 
field, including LkHa 198. Such spectra can indicate 
whether line emission is optically thick or self-absorbed 
(from the line shape) or if high- velocity gas is associated 
with a given position (though the presence of emission in 
the line wings). 

Figure [H presents the 12 CO J = 1 — high spatial and 
spectral resolution data from the combined FCRAO and 
BIMA array data cubes toward each of the four sources 
in the LkHa 198 region. The spectral resolution of 0.127 
km s _1 is a modest f actor of two improvement from the 
work of lCanto et al.l (|l984l ). Strong self-absorption is ev- 
ident toward each of the HAeBe objects of Figured! in- 
dicating that each remains embedded within the nascent 
molecular cloud. Line wings are detected toward both 
LkHa 198 and LkHa 198-mm, while there is a sugges- 
tion of some red-shifted higher velocity emission toward 
LkHa 198-IR, although this feature could be associated 
with the ambient cloud material since it arises in the 
combined spectra, but is not apparent in BIMA array 
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Fig. 4. — Spectra toward individual sources in the LkHa 198 
region from FCRAO and FCRAO+BIMA array data. The spectral 
resolution is 0.127 km s _1 . Grey, red and blue spectra represent 
12 CO J = 1-0, 13 CO J = 1-0 and CS J = 2-1 observations from 
FCRAO which were regridded to the spatial and spectral scales of 
the FCRAO+BIMA ^CO J = 1 - map. The intensity scale 
at left is for the single dish 12 CO J = 1 - emission. The 13 CO 
J = 1 — and CS J = 2 — 1 emission have been multiplied by factors 
of two and twenty. The black spectra are taken from the combined 
12 CO J = 1 — map. The intensity scale is shown at right. Self- 
absorption is noted toward each source in the field in both the 
single dish and combined 12 CO J = 1 — spectra. The strongest 
evidence of high velocity line wings is observed toward LkHa 198 
and LkHai 198-mm. Some red-shifted emission is observed toward 
LkHa 198-IR. 

data alone. 

Our FCRAO spectral line data enables us to assess 
whether tracers of dense (N(H2) ~ 10 5 cm -3 ) gas such 
as CS J = 2 — 1 and N2H + J = 1 — are present in the 
environment around LkHa 198. Figure [H shows 12 CO 
J = 1 - 0, 13 CO J = 1 - and CS J = 2 - 1 emission 
from FCRAO (integrated from a regridded map to match 
the BIMA array pixel scale, map size and velocity reso- 
lution). N2H + J = 1 — was not detected anywhere in 
the observed field. The CS J — 2 — 1 emission is peaked 
toward the self-absorption of the spectra derived from 
the high spatial resolution combined data. This is not 
surprising, since the CS J = 2 — 1, despite the limited 
spatial resolution of the maps, indicates the presence of 
relatively high density material (~ 10 4 cm -3 ). For the 
two northern sources, the CS J = 2 — 1 peak also aligns 
well with the peak of the 13 CO J = 1 — and the ab- 
sorption feature of single dish 12 CO J = 1 — 0, at 0.3 km 
s -1 . For the two southern sources, the 13 CO J = 1 — 
peak appears somewhat blue-shifted compared to the CS 
J = 2 — 1 peak, with the greatest offset toward LkHa 198. 

Ideally, we could use the high spatial resolution to solve 
for the V(LSR) as a function of position across the re- 
gion. However, the 12 CO J = 1 - emission spectra are 
too complex to be easily fit for single V(LSR) solutions. 
The CS J = 2 — 1 emission would be ideal, but for this 
species we have only FCRAO data. These data show no 
significant variation in peak velocity position across the 
FCRAO maps (i.e., 4'). 

A comparison of the 12 CO J = 1 — spectra from the 
single dish alone and the combined map (notably with 
different beams) indicates that the line wings are much 
more prominent in the spectra in which the BIMA ar- 
ray data are included. High resolution data enable us to 



discriminate between sources, e.g., clarifying the absence 
of line wings toward LkHa 198-IR. In addition, these 
data have increased sensitivity to outflow features which 
are significantly beam-diluted in the single dish maps. 
LkHa 198 has extensive line wings shifted blueward and 
redward of the source's velocity. Toward LkHa 198-mm 
and V376 Cas, the red-shifted emission is more promi- 
nent than the blue-shifted emission in the combined data 
but the opposite appears to be true in the single dish 
data, perhaps indicating that the blue lobe of the out- 
flow is more diffuse than its red counterpart. We also 
note that evidence for a blue wing is apparent in the 
13 CO J = 1 - spectrum toward V376 Cas. 

3.1.3. 2.6 mm Continuum 

No 2.6 mm continuum emission was detected toward 
the LkHa 198 region. This is consistent with the ab- 
sence of a detection at 2.7 mm of dust emission towar d 
LkHa 198 and V376 Cas bv iDi Francesco et all (|1997f ). 
Even in regions of 12 CO J = 1 - and N 2 H+ J = 1 - 
emission, no continuum emission is evident down to a 
threshold of 57 mJy bea m -1 , which is our 3 a limit. 
IDi Francesco et all (|1997l ) did detect 2.7 mm emission 
toward LkHa 198-mm using the Plateau de Bure IRAM 
interferometer. They detected a peak flux density of 
S v = 13 mJy beam -1 with a resolution of 9.7" x 5.1" 
and an integrated flux density of 24 ± 6.9 mJy. Our rms 
level of 18.5 mJy beam -1 is insufficient to detect this flux 
density. Our lack of detections also indicates that there 
is no significant free-free emission in the region. 

3.2. Individual Objects 

We summarize some general findings toward each 
source in our two fields in Table [H Some sources are as- 
sociated with CO emission, but not definitively with an 
outflow. All four sources in the LkHa 198 field are asso- 
ciated with local 12 CO J = 1 — emission, as revealed 
by the strong self absorption in each of the spectra of 
Figure ffl 13 CO J = 1 - and CS J = 2 - 1 are also 
detected toward each of the sources. 

3.2.1. LkHa 198 

Strong evidence for a blue-shifted cavity centered on 
LkHa 198 at a position angle 5 of ~ 120° is evident in the 
integrated channel maps of Figure O The same emission 
cavity is also obvious in the BIMA array data of Fig- 
ure 03 This emission is likely associate d with the cav- 
ity traced by near-IR I-band emission ([Corcoran et al.l 
119951 : lLagage et all 119931 : iLeinert et all [1991|). Theze- 
roth moment map of Figure [T] also highlights a weaker 
blue-shifted feature to the southwest of the LkHa 198 
sources (at a position angle of ~ 250° E of N). This 
could be a portion of a second outflow centered on LkHa 
198. 

The appearance of a single extensive blue lobe to the 
northwest and driven by LkHa 198 is show n to be a result 
of res olution limitations in earlier data. Nak ano et al.l 
(|199Q ) indicated that the velocity width of the blue lobe 
was significantly larger than the width of the red lobe 
they detected. Their single blue lobe is a combination 

5 We define the position angle according to the component of 
the outflow which is most clearly defined. 
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TABLE 4 
LkHo; 198 & Surrounding Objects 



Source 


Other ID 


RA 


DEC 


Group a 


Associated 


# Associated 


PA 






(J2000) 


(J2000) 




CO 


Outflows 


(° E of N) 


LkHa 198 


V633 Cas 


00:11:25.97 


+58:49:29.1 


II 


yes 


2 


120°, 70° 


LkHa 198-IR 


LkHa 198 B 


00:11:26.22 


+58:49:34.58 




yes 


none 




LkHa 198-mm 




00:11:24.51 


+58:49:42.81 




yes 


1 


350° 


V376 Cas 




00:11:26.21 


+58:50:04.1 


II 


yes 


1 


70° 



Note. — a Classification of IHillenbrand et al.l JT992|). 

of the blue-shifted emission driven from V376 Cas and 
the young, embedded source LkHa 198- mm. However, 
the d ouble-peaked, red-shifted emission of iNakano et al.l 
(|l99Qh persists even in our higher spatial resolution data. 
The peak observed to the north of LkHa 198 may arise 
from red-shifted emission from a LkHa 198-driven out- 
flow in combination with emission from an outflow driven 
by LkHa 198-mm, while the red-shifted emission to the 
south of LkHa 198 indeed arises from an outflow driven 
by that source. We suggest that at least one outflow 
is being driven from LkHa 198, based on the extended 
blue-shifted cavity. LkHa 198-IR has b een identified as 
the so urce associated with the cavity by ICorcoran et al.l 
(1995), but our CO J = 1 — observations suggest that 
LkHa 198 is driving at least one outflow to the southeast. 

Figure [1] compares the high resolution 12 CO J = 1 — 
da ta with the polarize d intensity at H-band measured 
by lPerrin et al.l (|2QQ4al ). The bipolar nebulosity centered 
on LkHa 198 in H-band emission has no counterpart in 
12 CO J = 1 - from BIMA array data. In fact, the 
high resolution components of the red- and blue-shifted 
12 CO J — 1 — emission is oriented primarily east and 
west of LkHa 198. The polarized emission should locate 
the circumstellar dust around the star; in this case, the 
north-south alignment could indicate the region of the 
star where the envelope has not yet been eliminated by 
the action of the outflow(s). 

3.2.2. LkHa 198-IR 

Based on Figure[TJ there is little evidence for any signif- 
icant outflow centered on LkHa 198-IR. Figure [H shows 
that there is no high velocity blue wing associated with 
the infrared source and only a suggestion of slight red- 
shifted emission. However, like LkHa 198, LkHa 198-IR 
shows strong self-absorption of the CO line, meaning that 
circumstellar gas is still associated with this object. 

The "cometary" structure detected in intensity and po- 
larization in the near-IR has no detected counterpart in 
12 CO J = 1 — emission. It is difficult to reco ncile the 
HH objects A and B from ICorcoran et ail (|l995f ) as orig- 
inating from LkHa 198-IR without any hint of a 12 CO 
J = 1 — outflow in our data. In fact, there is a hint of 
a red-shifted line wing in the 12 CO J = 1 — spectrum 
toward LkHa 198-IR (see Figure Hj). This could indicate 
a large-scale outflow contribution or ambient material 
contributions to the spectrum. 

3.2.3. LkHa 198-mm 

A roughly north-south outflow (with blue-shifted emis- 
sion at a position angle of ~ 350°) from LkHa 198-mm is 
strongly suggested by Figures [T] and El The blue-shifted 
component is particularly apparent in Figure [2] and ap- 



pears quite localized to the source. This is also suggested 
by the BIMA array data alone in Figure [3] although we 
note that there is significant emission missing from this 
map. The red lobes of the outflows of LkHa 198 and 
LkHa 198-mm clearly overlap. It is certain that part 
of the 12 CO J = 1 — red-shifted emission detected by 
ICanto et al.l (|l984f ) originates from this source and not 
LkHa 198. 

Figured! shows the self- absorption present in all sources 
in the region, indicating that substantial circumstellar 
gas remains. Figure [4] also shows that the red lobe from 
LkHa 198-mm is much stronger than the blue lobe. The 
fact that the outflow emission appears to lie so close to 
the driving source may indicate either that the outflow is 
almost entirely along the line of sight or that the driving 
source is very young. The absence of very high veloc- 
ity line wings argues against a line-of-sight alignment. 
Superior resolution will be required to disentangle the 
outflow of LkHa 198-mm from the other outflows in the 
region and the ambient cloud. We note that the ori- 
entation of the outflow is completely consistent with a 
series of HH objects observed a t position angles of 160° 
and 340° (C orcoran et al.lll995h which were attributed to 
LkHa 198, but are in fact more symmetric about LkHa 
198-mm. 

3.2.4. V376 Cas 

IHillenbrand et all 0992) categorize this source as a 
Group II HAeBe star, meaning that its SED is flat or ris- 
ing toward longer wavelengths. We therefore expect that 
there is a substantial amount of cool material around 
the star which reprocesses more of the total radiation 
than could be explained by a flat disk. It is therefore 
not surprising that the 12 CO J = 1 — spectrum of this 
source shows evidence of self- absorption in Figured! Ev- 
idence for nebulosity has been prev iously noted to the 
west of V376 Cas b y several authors ([Leinert et al.lll99ll: 
iPiirola et all 119921 : ISmith et al.l 12004^ iLeinert et all 
(|1991[ ) used polarized imaging techniques to detect the 
elongation of this nebulosity (120°) and clearly identify 
V376 Cas as the illumination source of the nebula. 

Figure [T] demonstrates the presence of extended, over- 
lapping blue- and red-shifted emission around this source 
at a position angle of ~ 70°. The BIMA array data alone, 
shown in Figure 02 suggest that the strong blue-shifted 
emission peak noted northwest of LkHa 198 in single- 
dish maps of the region may be dominated by emission 
from V376 Cas. This outflow orientation seems consis- 
tent with the nebulosity observed in the optical, given 
the difference scales to which those data and our molec- 
ular lines maps sample. The blue outflow is also coinci- 
dent with the known HH objects to its east, designated 
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Fig. 5. — Left: The greyscale is the zeroth moment over the entire range of CO emission (0 to 18 km s _1 ) of the combined FCRAO and 
BIMA array data toward the region centered on LkHa 225S (in units of Jy beam -1 km s -1 ). Optically visible HAeBe stars are marked 
by stars, while the embedded submillimeter source is indicated by a square. Contours trace blue emission from to 4 km s _1 and red 
emission from 12 to 18 km s _1 . The blue contour levels range from 4.5 to 14.4 Jy beam _1 km s _1 in steps of 0.9 Jy beam -1 km s _1 . The 
red contour levels range from 6.3 to 24.3 in steps of 1.8 Jy beam" 1 km s _1 . The resolution is 5.7" x 4.0", or 5700 x 4000 AU at 1 kpc. 
Two distinct blue-shifted features, "A" and "B" are observed. Position-velocity slices are taken through both features to align with the 
LkHa 225S (solid lines) and LkHo; 225N (dashed lines) sources. The box shows the field of view of the uniform-weighted image at right. 
Right: A combined map generated with the uniformly- weighted BIMA array data yields a resolution of 2.87" x 2.29". The symbols and 
velocity ranges of the zeroth moment maps are the same as for the naturally-weighted image. The blue contour levels range from 5 to 10 
Jy beam -1 km s _1 in steps of 1 Jy beam -1 km s -1 , and the red contour levels ra nge from 10 to 20 Jy beam -1 km s -1 in steps of 2 Jy 
beam -1 km s -1 . For comparison, the polarized intensity in H-band as measured by Perrin et al. (2004b) is shown in the inset. The LkHa 
225S and LkHa 225N sources are marked on this image. 



HH162 (see Fig. 1 of [Corcoran et al.lll995f ). The dom- 
inance of the blue-shifted emission to the east and the 
lack of comparable red-shifted emission to the west is 
entirely consistent with V376 Cas being located at the 
eastern edge of a small dark cloud core (see Fig. 1 in 
iSmith et al.l[2QQ3 ). 

4. LkRa 225S AND ENVIRONS 

Another example of an HAeBe star in a group is LkHa 
225. LkHa 225 comprises two sou rces separated by 
approximately 5" (|Aspin et al.lll994ft : LkHa 225N and 
LkHa 225S (V* V1318 Cyg N and V* V 1318 Cyg S) 
at a distance of 1 kpc (| Canto et al.l fl984 after iHerbid 
Il960f ). At wavelengths longward of 5 /mi, LkHa 225 
is the dominant source in th e region (L ~ I6OOL0, 
iThe. de Winter fc Pere3ll994f ). A third peak observed 
by lAspin et al.l (|1994| ) in the optical arises due to neb- 
ulosity between the two sources. This peak is coinci- 
dent with the comma shaped reflection neb ula extending 
from Lk H a 225N to L k Ra 22 5S detected bv lPerrin et~al] 
(|2004bD . iPalla et all (|1995[ ) detect two sources in the 
K-band with no e vidence of the third source noted by 
lAspin et al.l (|l994f ) and also detect an H2O maser co- 
incident with LkHa 225S lying at the center of a CO 
outflow. The detection of continu um peaks at the po- 
sition of LkRa 225S at 1.3 mm (|Henning et all I1998D 
and in high resolution imag i ng at 2.7 mm an d 3.1 
mm flPi Francesco et al.lll997t lLoonev et al.l 12006) con- 



firm that LkHa 225S is deeply embedded and likely 
the youngest star in the region. LkHa 225N is likely 
a pre-main sequen ce star. Near-IR stellar spectra 
([Davies et al.ll200T| ) show no absorption lines, only emis- 
sion lines, making the identification of the spectral types 
of both stars difficult. Two other well-studied HAeBe 
stars are in the same region: BD +40 4124 (V1685 Cyg, 
MWC 340, IRAS 20187+4111) and LkHa 224 (V1686 
Cyg), located approximately 50" to the northwest and 
15" to the wes t , resp ectively, from the LkHa 225 sources. 
lLoonev et all (|2006f ) show that the LkHa 225 sources 
lie in the brightest molecular cloud core (as traced by 
BIMA array observations of CS J = 2 — 1) in an elon- 
gated molecular filament. CS J = 2 — 1 emission at 2.5" 
resolution is strongly peaked on LkHa 225S. 



4.1. 



4.1.1. 



Observational Results 

Moment and Channel Maps 

Figure [5] shows the zeroth moment map summed from 
to 15 km s _1 , and the distribution of emission blue- 
ward and r edward of the sou rce's LSR systemic velocity, 
7.9 km s- 1 (jPalla et al.ll!995h . The 12 CO J = 1-0 emis- 
sion peaks around LkHa 225S. This concentration is un- 
surprising since a submillimeter peak is also detected at 
its position, indicating that the associated source is still 
deeply embedded and very young. The uniform- weighted 
image shows clearly that the CO peaks do not coincide 
with the positions of the potential driving sources, in- 
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Fig. 6. — Binned velocity channels of 12 CO J = 1 — toward the LkHo; 225S region are shown in greyscale and contours over 2.03 km 
s — 1 bins. The mean velocity per channel is noted in km s — 1 . The greyscale range is common to all channels and is shown at right (Jy 
beam -1 ). The contours are plotted at 20 and 40 to 360 a in steps of 40 a where a = 0.45 Jy beam -1 is the rms level in a single channel 
before binning. The sources are marked by asterisks and a square as in Figure GD 

4Q m 

dicating that the 12 CO J = 1 — emission associated 
with these sources is primarily in outflowing gas. Fig- 
ure [6] shows channel maps of 12 CO J = 1 — emission. 
The 12 CO J = 1 — morphology is very complex as ev- 
idenced by the variation of structure with velocity. The 
most intriguing structure is an apparent ring of emission, 
centred at a J20 oo = 20 h 20 m 30*0, S J2 ooo = +41°21 / 28 / /0, 
which persists over a large range of velocity (10 to 18 km 
s _1 ). A SIMBAD 6 search lists no objects at the center 
of the ring; the closest objects are the HAeBe stars in 
the field. An H2O maser is reported at the position of 
LkHa 225S, identified as the driving s ource of the 12 CO 
J = 2-1 outflow seen bv lPalla et all (|l995h at 15" reso- 
lution . The blue-shifted emission detected bv lPalla et al.l 
(|l995f ) lies just to the east of the peak of the red-shifted 
emission, implying that the outflow lies almost along the 
line-of-sight. Peaks on the ring are coincident with the 
positions of LkHa 224 and LkHa 225S, but LkHa 225N 
lies between 12 CO J = 1 — peaks. We surmise that 
the ring is the supposition of emission from several out- 
flows, dominated by the pole-on outflow, which could 
have a substantial opening angle and could appear as a 
cylindrical cavity. No obvious counterpart to the ring is 
observed in blue-shifted emission, but strong blue-shifted 
emission is concentrated at the position of LkHa 225S. 

To ascertain the physical origin of the "ring", inspec- 
tion of position- velocity (P — V) diagrams is critical. Fig- 
ure [71 shows the difference between two P — V slices, one 
taken centered on LkHa 225S and one taken at a position 
offset 24" to the south, with the same orientation. Taking 




6 Guest User, Canadian Astronomy Data Centre, which is oper- 
ated by the Dominion Astrophysical Observatory for the National 
Research Council of Canada's Herzberg Institute of Astrophysics. 
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Velocity (km/s) 

Fig. 7. — The top panel shows a position velocity diagram taken 
centred on LkHo; 225S sliced through the blue-shifted feature "A" 
of Fig. [5] The evidence for multiple clouds along the line of sight 
is indicated by the bands of emission at all angular positions. The 
bottom panel shows the same position velocity slice, but with a 
second slice, positioned off the source 24" to the southeast, sub- 
tracted from the on-source slice. The bulk of the cloud emission 
is removed, especially the fainter clouds. In this image, the most 
prominent feature is the emission extending from about km s — 1 
to 20 km s _1 at zero angular offset. We interpret this as evidence 
for a pole-on outflow, most likely centered on LkHa 225S. Blue- 
shifted emission is also detected to the northeast of LkHo; 225S. 

the difference of these two P — V diagrams crudely re- 
moves the large scale emission from three separate molec- 
ular cloud components along the line of sight in Cygnus 
(at roughly 3 km s _1 , 8 km s _1 and 13 km s _1 ) and 
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Fig. 8. — Position-velocity diagrams taken through the two blue features marked A and B in Figure[5]are shown aligned through both the 
LkHa: 225S and LkHa: 225N sources. It is not possible to discriminate whether feature "A" arises from LkHa 225S or LkHa 225N. Feature 
B shows virtually the same profile through both sources. Since this feature is aligned with the offset between the sources, discrimination 
is naturally difficult. There is still some evidence for the pole-on outflow of LkHa 225S (see Fig. in slices through LkHa 225N. 



clearly reveals the presence of emission at the 0" offset 
across both red and blue-shifted velocities (from to 20 
km s _1 ) with respect to the source. The P — V diagram 
produces this feature no matter what the orientation of 
the slice centered on LkHa 225S. We interpret this emis- 
sion as arising from an outflow driven from LkHa 225S 
but observed almost pole-on. This i s consistent with the 
12 CO J = 2-1 emission detected bv lPalla et all (|1995j ) in 
which they find strongly overlapping blue and red lobes 
in the single detected outflow. The blue-shifted emission 
is offset slightly to the east compared to the red lobe, 
which is also consistent with the plane-of-sky emission 
shown in Figure but i s in the opposite sense to the 
outflow lobes observed bv lPalla et al.l (|l995l ). 

Figures [5] and [6] provide further clarification in identify- 
ing the driving sources of the outflows with components 
in the plane of the sky. In addition to the pole-on out- 
flow, a second outflow appears well-centered on LkHa 
225, based on the blue lobe of emission detected south- 
east of the two sources at a position angle of ~ 145° (E 
of N). Another blue peak is detected to the northeast 
of the cluster of sources (at ~ 45°), suggesting a third 
outflow is present in the system. There is no evidence of 
coherence in the blue-shifted emission peaks to suggest 
they are part of a larger, blue-shifted counterpart to the 
red-shifted "ring", i.e., part of the same outflow. 

Figure [8] shows the P — V diagrams across two po- 
sition angles suggested by the orientation of the blue- 
shifted emission detected in the zeroth moment and chan- 
nel maps (labeled "A" and "B" in Fig. Ej). Along these 
slices, there is evidence for red-shifted components com- 
plementary to the blue-shifted emission. The red-shifted 
emission from these two potential outflows appears jux- 



taposed on the zeroth moment map of Figure El lying 
coincident with the red-shifted gas from the pole-on out- 
flow and creating a ring structure which persists over a 
range of velocities, although this emission appears signif- 
icantly closer to the driving sources than the blue lobes. 

4.1.2. Spectra 

Using FCRAO, we have detected 12 CO J = 1-0, 13 CO 
J = 1 - 0, CS J = 2 - 1 and N 2 H+ J = 1 - emission 
from the LkHa 225S region. The 12 CO J = 1 - 0, 13 CO 
J = 1 — and CS J = 2 — 1 spectra all indicate multiple 
emission peaks along the line of sight and are shown to- 
ward each of the four sources in Figure [9l We also show 
the combined 12 CO J = 1 — spectra at high spatial 
and spectral resolution. We show an N 2 H + J = 1 — 
spectrum toward LkHa 225S in Figure [TUl Strong evi- 
dence for outflowing gas is detected toward LkHa 225S 
and LkHa 225N, but none is detected toward BD +40 
4124 and only evidence for red-shifted emission is seen at 
the position of LkHa 224. We interpret the high veloc- 
ity emission toward LkHa 224 as arising from the out- 
flows generated by the LkHa 225 sources. Less ambient 
CO appears to be present in this region compared to 
the environment around LkHa 198, since no significant 
self-absorption is observed toward any of the sources in 
Figure [9l Interpretation of these spectra is further com- 
plicated by the presence of multiple clouds along the line 
of sight. Peaks are observed at five different velocities 
in the spectra (1.6, 3.2, 7.9, 12.8 and 14.5 km s _1 as 
identified by eye in Figure [9|). The primary cloud, at 7.9 
km s _1 , appears to contain all the HAeBe stars. A sec- 
ond cloud, at 12.8 km s _1 , appears to contain significant 
amounts of 12 CO J = 1 - and 13 CO J = 1 - 0, but no 
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Fig. 9. — Spectra toward individual sources in the LkHa 225S 
region from FCRAO and FCRAO+BIMA array data. The spectral 
resolution is 0.254 km s _1 . Grey, red and blue spectra represent 
12CO J = 1-0, 13 CO J = l- 0andCSJ = 2- l observations 
from FCRAO which were regridded to the spatial and spectral 
scales of the FCRAO+BIMA 12 CO J = 1 - map. The intensity 
scale at left is for the single dish 12 CO J = 1 — emission. The 
13 CO J = 1 — and CS J = 2 — 1 emission have been multiplied 
by factors of two and twenty. The black spectra are taken from 
the combined 12 CO J = 1 - map. The intensity scale is shown 
at right. Unlike the sources in the LkHai 198 region, those in 
Cygnus do not exhibit significant self-absorption. LkHai 225S and 
LkHa: 225N show the most extensive line wings. BD +40 4124 
does not show any evidence of outflow while LkHa 224 shows no 
high velocity blue-shifted emission, but extensive red-shifted line 
wings which could be associated with another source in the field. 
The lines below the plots indicate 1.6, 3.2, 7.9, 12.8 and 14.5 km 
s _1 . These peaks each represent individual clouds along the line 
of sight. 
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Fig. 10. — Spectrum of N2H+ J = 1 — emission toward LkHa: 
225S, indicating the presence of dense gas. 

CS J = 2 — 1 was detected at this velocity. The evidence 
for clouds at the remaining velocities is based on 12 CO 
J = 1 - only. 

The resolution of FCRAO is such that we cannot be 
sure that all four sources are associated with dense gas. 
As the most embedded object, LkHa 225S is the most 
likely source to be associated with dense gas. How- 
ever, BD +40 4124, which is the furthest removed from 
LkHa 225S, lies only about 45" distant. This is only one 
FCRAO beamwidth at the frequency of the 13 CO J = 



1 — emission, and substantially less than a beamwidth in 
CS J = 2 - 1 and N 2 H+ J = 1 - 0. Higher resolution ob- 
servations can isolate the concentrations of these species 
and enable kinematic study of the associated cores. 

4.1.3. 2.6 mm Continuum 

No signficiant continuum was detected toward any of 
the sources in the field containing LkHa 225S. Our sensi- 
tivity is rather low, given the limited bandwidth available 
for continuum observations. We can put an upper limit 
of 51 mJy beam -1 on the continuum flux on compact 
sources in the region. Despite the strong CO and N 2 H + 
J = 1 — detections, we see no flux t o attribute to either 
dust continuum or free-free emission. Di Francesco et al. 
(1997) detected 2.7 mm emission from LkHa 225S using 
the Plateau de Bure interferometer with a resolution of 
6.8" x 4.6". They detected a peak flux density of 23 mJy 
beam -1 and a total flux density of 24 ± 6.7 mJy. This 
flux would not have been detectable at a 2 <j level with 
the sensitivity (a = 17 mJy beam -1 rms) in our map. 

4.2. Individual Sources 

Table [5] summarizes our findings regarding 12 CO J = 
1 — and outflows in this region. Self-absorption is not 
detected toward any sources in this field. However, we 
do detect tracers of dense gas, including 13 CO J = 1 — 0, 
CS J = 2 - 1 and N 2 H+ J = 1 - in this cloud. We note 
also that our interpretations of outflows from LkHa 225S 
and LkHa 225N are still resolution-limited. Therefore, 
we discuss these sources together below. 

4.2.1. LkHa 225S and LkHa 225N 

LkHa 225S is found to have evidence for up to three 
associated outflows. The most definitive connection is 
based on evidence for a nearly pole-on outflow clearly 
seen in the P — V diagram of Figure [3 Evidence for two 
additional outflows is seen in Figure [5] in the form of two 
blue-shifted emission features. The solid lines show the 
axes formed by the features, "A" and "B", with LkHa 
225S and both axes align with red-shifted features seen 
on the other side of the potential driving source. How- 
ever, LkHa 225S is separated from LkHa 225N by only 
5", which is comparable to the major axis of our synthe- 
sized beam (naturally- weighted). This makes it difficult 
to associate decisively either of these two plane-of-sky 
outflows with LkHa 225S. Based on alignment, a strong 
case can be made that LkHa 225N is driving one or both 
of these outflows (see the dashed lines of Figured]). 

We assert that the apparent ring of Figure [5] has coher- 
ent structure, as it is dominated by the red-shifted emis- 
sion of the pole-on outflow with contributions from the 
other outflows originating from the LkHa 225 sources, 
one oriented at approximately 45-55° E of N, the second 
at about 145-160° E of N and the third almost pole-on to- 
ward us. No evidence of the ring morphology is observed 
in pre vious detections o f high- velocity gas in the region 
(e.g., iPalla et "all 1 19951 ) with resolution only a factor of 
three poorer than our data. 

As we noted above, the interpretation of the ring as 
the red-shifted counterparts to the more definitive blue 
lobes on the plane of the sky suggests that the red lobes 
have no t expanded to the sam e spatial extent. The recent 
work of lLoonev et al.l (|2006l ) suggests a potential expla- 
nation. Their high resolution BIMA array observations 
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TABLE 5 

LkHq; 225S & SURROUNDING OBJECTS 



Source Other ID RA DEC Group a Associated # Associated PA 

(J2000) (J2000) CO Outflows (° E of N) 

LkHa 225S V1318S Cyg 20:20:30.6 +41:21:25.5 ... yes l-3 b l.o.s. c , 45°, 145° 

LkHa 225N V1318N Cyg 20:20:30.6 +41:21:30.5 ... maybe 0-2 b 55°, 160° 

BD +40 4124 V1685 Cyg 20:20:28.25 +41:21:51.6 I yes none 

LkHa 224 V1686 Cyg 20:20:29.26 +41:21:28.6 I yes none 



Note. — a Classification of Hillenbrand et al. (1992). b The driving source of two outflows cannot be uniquely 
identified as LkHa 225S or LkHa 225N due to remaining resolution limitations. c line-of-sight. The outflow is almost 
pole-on, with little component in the plane of the sky. 



of CS J = 2 — 1 show that the CS J = 2 — 1 traces out a 
dense filament of emission to the northwest of the LkHa 
225 sources. The presence of denser gas in this direction 
may be impeding the red-shifted component of the out- 
flow from extending into the cloud along these position 
angles. 

4.2.2. BD +40 4124 

BD +40 4124 is a Group I HAeBe star as classified 
bv [Hillenbra nd et all (11992 ). This means that there is 
no remaining circumstellar emission which cannot be ac- 
counted for by a disk. 12 CO J = 1 - 0, 13 CO J = 1 - 
and CS J = 2 - 1 are all detected toward BD +40 4124 
based on the spectra of Figure [U This emission is likely 
associated with the ambient molecular clouds along the 
line of sight to Cygnus. The lack of self-absorption indi- 
cates that relatively little gas remains around this source, 
and there is no significant evidence of line wings in the 
spectra. 

4.2.3. LkHa 224 

LkHq 224 is also classi fied as a Group I HAeBe star by 
iHillenbrand et all (|1992j ). The spectrum of LkHa 224 is 
dominated by the multiple peaks associated with cloud 
emission along the line of sight. The combined 12 CO 
J = 1 — spectrum does show evidence of a red-shifted 
line wing however. No outflow emission appears focussed 
on this source based on the channel maps of Figure^ and 
no blue-shifted emission is in the vicinity of this source in 
the zeroth moment map of Figure [U Given its position 
relative to the blue-shifted emission peaks, we surmise 
that the red-shifted emission observed at the position of 
this source is part of one of the outflows originating from 
the LkHa 225 sources. 

5. CONCLUSIONS 

We find that the single bipolar outflow detected in the 
LkHa 198 region in single-dish maps is a combination 
of the red-shifted emission of the LkHa 198 and LkHa 
198-mm outflows and the blue-shifted components of the 
outflows from LkHa 198-mm and V376 Cas. In all, we 
find evidence for four outflows in the LkHa 198 region 
using high-resolution 12 CO J = 1 — observations. Ev- 
idence for a cavity is observed in the LkHa 198 system 
along the same axis as the elliptical cavity detected in 
optical and IR emission. We see no evidence for circum- 
stellar material based on 12 CO J = 1 — emission as- 
sociated with LkHa 198-IR, contrary to the results from 
near-IR pola rimetry which ass ociate dust emission with 
that source (|Perrin et al.ll2004al ) and observations of HH 



objects, which identifi ed LkHa 198-IR as a driving source 
([Corcoran et alJll995[ ). 

Our ability to discriminate between detected 12 CO 
J — 1 — outflows and measure independent masses and 
momentum is limited by the single pointing of the BIMA 
array maps (Figure [3] clearly shows that we have not 
mapped the outflows in their entirety) and spatial reso- 
lution. The outflow emission is still entangled and the 
cloud emission is also a contaminant. Ideally, mosaicked 
observations with various array configurations (includ- 
ing very close baselines) will be obtained to maximize 
the detection of the outflow emission, but minimize the 
detection of emission from the ambient cloud. 

The 12 CO J = 1 - emission from LkHa 225S is 
equally complex. In all, we find evidence for up to three 
outflows in this region, all originating from the LkHa 225 
sources. Large-scale cloud emission dominates the field, 
and multiple clouds along the line of sight contaminate 
the spectra. One is definitively associated with LkHa 
225S and is almost entirely along the line of sight. The 
other outflows lie at angles of 35° and 135° east of north 
projected onto the plane of the sky. We suggest that 
the ring of red-shifted emission is in fact a juxtaposi- 
tion of the red-shifted lobes of these two outflows, where 
their expansion has been hindered due to a dense molec- 
ular cloud lying to the northwest, with the cavity of the 
pole-on outflow (a coherent structure). As in the LkHa 
198 region, higher resolution, mosaicked maps would help 
disentangle the outflows from the LkHa 225 sources. 

Estimation of the mass in these outflows is not easy 
based on the data presented here because they can only 
be resolved within the BIMA array antenna field of view, 
whereas the outflows certainly extend to larger spatial 
scales. In several cases, emission from several outflows is 
clearly combined, further complicating the independent 
assessments of masses. 

Finally, we find evidence of higher density gas in both 
regions, with 13 CO J = 1 - and CS J = 2 - 1 detected 
toward all eight sources and N2H + J = 1 — detected 
in the objects in Cygnus. Localizing the precise sources 
of this emission requires higher resolution observations, 
since the separation of sources is comparable or smaller 
than the FCRAO beam. 
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